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ABSTRACT 

We explore the dynamics of small planetesimals coexisting with massive protoplanetary cores in a 
gaseous nebula. Gas drag strongly affects the motion of small bodies leading to the decay of their 
eccentricities and inclinations, which are excited by the gravity of protoplanetary cores. Drag acting 
on larger (> 1 km), high velocity planetesimals causes a mere reduction of their average random 
velocity. By contrast, drag qualitatively changes the dynamics of smaller (< 0.1 — 1 km), low velocity 
objects: (1) small planetesimals sediment towards the midplanc of the nebula forming vertically thin 
subdisk; (2) their random velocities rapidly decay between successive passages of the cores and, as a 
result, encounters with cores typically occur at the minimum relative velocity allowed by the shear 
in the disk. This leads to a drastic increase in the accretion rate of small planetesimals by the 
protoplanetary cores, allowing cores to grow faster than expected in the simple oligarchic picture, 
provided that the population of small planetesimals contains more than roughly 1% of the solid mass 
in the nebula. Fragmentation of larger planetesimals (> 1 km) in energetic collisions triggered by 
the gravitational scattering by cores can easily channel this amount of material into small bodies on 
reasonable timescales (< 1 Myr in the outer Solar System), providing a means for the rapid growth 
(within several Myr at 30 AU) of rather massive protoplanetary cores. Effects of inelastic collisions 
between planetesimals and presence of multiple protoplanetary cores are discussed. 
Subject headings: planetary systems: formation — solar system: formation — Kuiper Belt 



1. INTRODUCTION. 

Formation of terrestrial planets and solid cores of gi- 
ant planets is thought to proceed via the gravity-assisted 
merging of a large number of planetesimals — solid bod- 
ies with initial sizes of roughly several kilometers. De- 
spite a considerable progress made in this field since the 
pioneering works of Safronov (1969), a number of impor- 
tant problems still remain unsolved. One of the most se- 
rious questions has to do with the time needed for planets 
to complete their growth to present sizes. In the frame- 
work of conventional theory this time is rather long, es- 
pecially in the outer parts of the protoplanetary nebula 
Qt 10 8 — 10 9 y r )j an d it is likely that the gaseous compo- 
nent of the nebula dissipates much earlier (in < 10 6 — 10 7 
yr) . This would make it very hard for the giant planets in 
our Solar System to accrete their huge gaseous envelopes 
via core instability (Mizuno 1980) which is otherwise con- 
sidered to be an attractive scenario. 

Wetherill & Stewart (1989) have identified a very rapid 
"runaway" regime of accretion of protoplanetary cores 
which at the time seemed like a solution of this problem. 
However, later on Ida & Makino (1993) and Kokubo & 
Ida (1996, 1998) have demonstrated that the runaway 
accretion would persist only through a rather limited in- 
terval of time and that the final growth of protoplan- 
etary cores to isolation (which corresponds to roughly 
10 26 - 10 27 g at 1 AU) would proceed in a slow manner, 
making the formation of cores of giant planets rather 
problematic. 

These studies usually implied that the gaseous com- 
ponent of the nebula plays only a secondary role in the 
planet formation process. Planetesimals were typically 
assumed to be rather massive (10 23 — 10 24 g) bodies 
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weakly affected by the gas (Kokubo & Ida 1998). This al- 
lows gravity to excite energetic random motions of plan- 
etesimals leading to diminishing the role of gravitational 
focusing and reduction of accretion efficiency. The pur- 
pose of this paper is to relax this assumption and to see 
what impact is incurred on the planet formation picture 
by allowing most of the planetesimals to be small (< 10 
km) bodies immersed in a gaseous environment. Such 
planetesimals would be appreciably affected by the gas 
drag and we will demonstrate that this can bring qual- 
itative changes to their dynamics in the vicinity of the 
protoplanetary cores and, consequently, to the behavior 
of the mass accretion rate of cores. 

Throughout this study we will use the following ap- 
proximation to the structure of the Minimum Mass Solar 
Nebula (MMSN): 
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where are the particulate and gas surface den- 

sities correspondingly, c s is the gas sound speed, and 
clau = a /(l AU) is a distance from the Sun a scaled 
by 1 AU. We will use the terms "protoplanetary core" 
and "protoplanetary embryo" interchangingly. Physical 
density of planetesimals p p is always assumed to be 1 g 
cm~ 3 . 

The paper is organized as follows: after a discussion 
of different gas drag regimes in |J2]we proceed to the de- 
scription of planetcsimal dynamics in the vicinity of pro- 
toplanetary cores in The inclination of small plan- 
etesimals, a question very important for this study, is 
explored in WS 3.41 The separation of different gas drag 
and dynamical regimes in different parts of the nebula 
is described in H3 3.51 and lower limits on the random ve- 
locities of planetesimals are obtained in 21 We dwell 
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upon the role of inelastic collisions between planetesi- 
mals in SjSl The role of small planetesimals in the growth 
of protoplanetary cores is studied in ^SJand some impor- 
tant consequences for the planet formation picture are 
discussed in SJ3 

2. SUMMARY OF DIFFERENT GAS DRAG REGIMES. 

Drag force acting on a body moving in a gaseous 
medium depends on the relative velocity of the body with 
respect to the gas v r and on the ratio of its radius r p to 
the molecular mean free path A (Whipple 1972; Weiden- 
schilling 1977). Whenever r p < A the Epstein drag law 
applies: 
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where f2 is the angular frequency of the disk, p g is the gas 
density, r p is the planetesimal size, and p p is the physical 
density of planetesimals (planetesimal random velocities 
are assumed to be subsonic). For the adopted MMSN 
model we estimate 
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for H2 collision cross section <jh 2 — 10~ 15 cm 2 (p is the 
mean molecular weight). Here n g is a number density 
of H2 molecules. It is important to notice that A in- 
creases very rapidly with the distance from the Sun, so 
that although only sub-cm particles can experience Ep- 
stein drag at 1 AU, at 30 AU this drag regime is valid 
even for rocks 100 m in size! 

A large spherical body with r p > A experiences a de- 
celeration of the form 
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Drag coefficient Cd is a function of the Reynolds num- 
ber Re = v r r p /v, where v is a kinematic viscosity. The 
viscosity we assume here is the molecular viscosity, i.e. 
v ~ Ac s /3; in the presence of anomalous sources of vis- 
cosity the expression for v has to be correspondingly ad- 
justed. 

For Re < 1 gas drag is in the Stokes regime (Landau & 
Lifshitz 1987), meaning that C D = 6vri?e _1 . For Re > 1 
drag coefficient becomes constant: Cjj ~ 0.7 (Weiden- 
schilling 1977). We neglect the more complicated behav- 
ior of C D for 1 < Re < 10 2 (Whipple 1972) and simply 
assume that Cd switches from one asymptotic behavior 
to the other at Reb — 20. Thus, we adopt 
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Separation of different gas drag regimes as a function of 
r p and v r is schematically presented in Figure Q 

Gas drag acting on a planetesimal moving on an elliptic 
and inclined orbit leads to the decay of planetesimal ec- 
centricity e, inclination i, and semimajor axis a (Adachi 
et al. 1976). For the calculation of the decay rate of the 
orbital elements it is important to take into account the 
sub-Keplerian angular velocity of the gas in the nebula 
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Fig. 1. — Separation of different gas drag regimes on r p / X - v r /c s 
plane. Dashed lines separate regions where different gas drag laws 
(indicated on a plot) apply. Planetesimals in the shaded part of the 
plot move supersonically with respect to the gas. Separation be- 
tween the quadratic and Stokes regimes (Re = -Re;,) is determined 
by v r /c s = (iieb/SXrp/A)- 1 . 



caused by the radial pressure support. This gives rise 

j to the relative 
One can easily demonstrate 



to additional azimuthal contribution Av g to the relative 



gas-planetesimal velocity, 
that (Whipple 1971) 
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For our adopted temperature profile Av g is independent 
of a. Whenever random velocity of planetesimal exceeds 
Av g this velocity offset provides only a small contribution 
to the relative gas-planetesimal velocity. In the opposite 
case relative velocity is dominated by Av g . 

Using ©-0 we can introduce a damping time td - 
typical time needed to decelerate planetesimals by the 
gas drag if their initial velocity with respect to the gas is 
Av g . Using equations ©, ©, & one obtains 
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for the quadratic, Stokes, and Epstein drag regimes re- 
spectively. 

When td <§C planetesimal motion is tightly coupled 
to that of the gas; critical planetesimal size r stop at which 
the transition from almost Keplerian motion to the sub- 
Keplerian gas rotation occurs can be determined from 
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the condition t d f« Using (JTJU-tEJ we find that 
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for the quadratic, Stokes, and Epstein drag regimes. 

Calculation of the orbital evolution of planetesimals 
bigger than r stop due to gas drag was performed by 
Adachi et al. (1976; see also Tanaka & Ida 1997) tak- 
ing into account the saturation of the relative gas- 
planetesimal velocity whenever the planetesimal random 
speed is < Av g . We quote their results for different 
drag regimes keeping only the most important contri- 
butions : 
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for the quadratic dependence Q of gas drag on velocity. 
The decay rates represented by equations flfi)) & fTTj) be- 
have differently for e, i < r/ and e, i > 77. In the first case 
the relative gas-planetesimal velocity is close to Av g ; as 
a result, orbital elements exponentially decay with time. 
In the second case v r is dominated by the planetesimal 
random motion and the decay time is inversely propor- 
tional to the random velocity v ss (e + i)Qa. 

For the linear dependence of gas drag on v r represented 
by the equations and © one finds (Adachi et al. 
1976) 
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where £ = 1 for r p 
(3/2)A/t> for A < r p 



< A (Epstein drag [3]), and £ 
i?e < i?e b (Stokes drag jSJ). 



3. SCATTERING BY MASSIVE PROTOPLANETARY CORES. 

Let us consider the scattering of a planetesimal by a 
single protoplanetary core moving on a circular and unin- 
clined orbit. The velocity perturbation which planetes- 
imal receives in the course of gravitational interaction 
depends on the separation of the planetesimal guiding 
center from the embryo's orbit and on the planetesimal 
approach velocity v a . 

It is well known (Henon & Petit 1986) that the Hill 
radius Rh of interacting bodies sets an important length 
scale for scattering in the disk. For two bodies with 
masses mi and m 2 Hill radius is defined as 
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1 We have set to unity all numerical factors which appear in the 
course of orbital averaging; however we kept initial factors since 
some of them (e.g. in 0) are significantly different from unity. 



In the case of scattering of planetesimals by a protoplan- 
etary embryo having mass M e Hill radius is determined 
solely by M e : R H = a(M e /M )V3. 

For our purposes the planetesimal approach velocity 
can be represented as v a ~ fl(ea + ia + Rh) and it is 
different from the horizontal random velocity = efla 
and vertical random velocity v z = ifla. Additional 
contribution to v a in the form of the shear across the 
Hill radius HRh (sometimes called the Hill velocity) ap- 
pears because of the differential rotation of the disk. 
Gravitational focusing of planetesimals by the embryo 
is determined by v a . Velocity v z determines the verti- 
cal thickness of planetesimal disk, thus it regulates the 
volume number density of planetesimals n p for a given 
surface mass density E p : n p = fl^ p /(v z m p ), where 
m p = (4-7T /i)p p Tp is a planetesimal mass. Vertical ve- 
locity is much smaller than either v a or the total random 
velocity v ~ + v z when i <C e. 

Whenever random velocity v is less than - 
the so-called shear-dominated regime (Ida 1990) — close 
approaches between gravitationally interacting bodies 
which might lead to their physical collision are only pos- 
sible for the orbital separation h ~ Rh- For h ^> Rh 
(distant encounters) orbital perturbations incurred in the 
course of scattering are small: changes of the random ve- 
locity Av and of the orbital separation Ah produced by 
a single scattering have magnitudes 

Av «fLR ff (R B /hf , AhK R H {R H /h) 2 . (21) 

These increments can be positive or negative depending 
on the orbital phases of bodies and their random veloci- 
ties before the encounter. The averages of these quanti- 
ties over epicyclic phases were calculated by Hasegawa & 
Nakazawa (1990) and we quote here their result for the 
eccentricity (it will be used in 
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In the case h ~ Rh, when interacting bodies enter 
their mutual Hill sphere, strong scattering takes place: 

Av w QR H , A(h 2 ) « R%, (23) 

with both increments being positive — orbits are repelled 
and random velocity ~ QRh is excited. At the same 
time, it is clear from purely geometrical considerations 
that for v <C flRn the change of the vertical velocity is 
proportional to the vertical projection of the random ve- 
locity perturbation, and is much smaller than the change 
in the horizontal velocity: 



Av z 
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(24) 



Since the total velocity increment Aw is of the order or 
bigger than the initial random velocity of planetesimals, 
scattering for h ~ Rh in the shear-dominated regime has 
a discrete character. 

Another possible dynamical state of planetesimals is 
the dispersion-dominated regime which takes place when- 
ever v > HRh- In this case the epicyclic excursions of 
planetesimals allow collisions between the bodies with or- 
bits separated by h < v/Cl; more distant bodies are again 
subject to only weak scattering. Moreover, even particles 
having close approaches experience on average only small 
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velocity perturbation compared to the pre-encounter ve- 
locity v a . This allows one to treat the scattering in 
the dispersion-dominated regime as a continuous process, 
unlike the discrete scattering in the shear-dominated 
regime. Another important feature of the dispersion- 
dominated case is that v z ~ Vh, or i ~ e, as a result 
of roughly three-dimensional nature of scattering in this 
regime which tends to isotropize highly anisotropic veloc- 
ity distributions (Ida, Kokubo, & Makino 1993; Rafikov 
2003c). Thus, it is enough to follow the evolution of only 
one component of planetesimal velocity (e.g. e) in the 
dispersion-dominated regime. 

3.1. Quadratic drag. 

Let's now proceed to considering planetesimal dynam- 
ics for different gas drag regimes. First, we explore the 
case of quadratic drag represented by equations , <|lfi|) , 
and i|17|l which is valid for r p > A, Re > i?e&. This situ- 
ation would most likely be realized in the inner parts of 
the protoplanetary nebula since A is rather small there 
making Epstein regime irrelevant and bringing Reynolds 
number to a high value. 

We first assume that scattering occurs in the 
dispersion-dominated regime and try to figure out under 
which conditions a steady state can be realized in this 
case. Gravitational scattering by the embryo increases 
planetesimal random energy at a rate 2 (Ida & Makino 
1993; Rafikov 2003b) 
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where In A ~ 1 is a Coulomb logarithm and A w 
(v/QRh) 3 for e - i > R H /a (Stewart & Ida 2000). Be- 
cause of the very weak dependence of In A on velocity we 
will set it to unity in our further discussion. 

Balancing this growth rate by the damping due to the 
gas drag (equation we find that 
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for v < Av g . (27) 

For further convenience we normalize the embryo's 
mass M e to a fiducial mass Mf and planetesimal radius 
to a fiducial radius r/ which are defined as 

M f = M r? 3 ~ 8 x 10 24 g a 3 j£, (28) 

r/ = .2^— = 0.2a^ ~ 150 m a"[/ 4 . (29) 
Pp c s p p 

These two parameters uniquely characterize the dynam- 
ics of the embryo-planetesimal scattering when the gas 

2 In the dispersion-dominated regime v a ~ V. 



drag dependence on the velocity is represented by equa- 
tion 10. Using this notation we conclude from Ij26(l that 
v > flRn and v > Av g if 

rp/r f > 1, (30) 

Me/M f > (r p /r f y 1/2 . (31) 

From H27|) we find that v > HRh and v < Av g if 

M e /M f > (r p /r f y 3 , (32) 

Me/M f < (r p /r f y 1/2 . (33) 

Thus, conditions l|30|l and l|32|) determine the region in 
the parameter space M e — r p in which planctcsimals are 
dispersion-dominated with respect to embryo and main- 
tain a steady-state velocity dispersion, see Figure 01 

Whenever both (|30() and (|32|l are not fulfilled, scat- 
tering of planetesimals proceeds in the shear-dominated 
regime. As we have previously mentioned, this type of 
interaction has a discrete nature — gas drag causes sig- 
nificant evolution of planetesimal velocity between close 
approaches to the embryo. We consider a planetesimal 
initially separated in semimajor axis from the embryo by 
h R H with initial random velocity v < flRn- Im- 
mediately after scattering v increases to ~ QRh, and 
the post-scattering value of planetesimal eccentricity is 

e ^R H /a = r l ( Me/Mf ) 1/3 - 

We will see in *I3 3.4l that the inclination of planetesi- 
mals interacting with the embryo in the shear-dominated 
regime is typically much smaller than their eccentricity. 
Then we can easily solve equation I jlfijl with initial con- 
dition e(0) = eo to obtain a general solution in the form 



e(t) « n 
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(34) 



For low-mass embryos, M e < Mf, planetesimal veloc- 
ity after shear-dominated scattering is always less than 
Av g , i.e eo < rj; as a result, we find from l|34() that v 
damps exponentially: 



e(t) i=a eo exp 
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Thus, eccentricity decays on a typical timescale td, see 
G3 & (El. 

For high-mass embryos, M e > Mf, post-encounter ve- 
locity is above Av g (eo > rj) and damping by the gas 
drag is very efficient: dv 2 /dt oc — v 3 , see equation {7J|. 
One can find from (|34() that in this case 

g(t) " (ry/eo) + nt V (r f /r p ) = ( V /e Q ) + (t/t d ) ' (36) 

for t < t c i- One finds that e drops to rj after t rs t^ 
independent of e . Beyond that point relative gas- 
planetesimal velocity is set by Av g , thus for t > td ec- 
centricity decays exponentially on a timescale td, anal- 
ogous to H35|) . It also follows from H3(jfl that signifi- 
cant reduction of eccentricity after scattering by the high 
mass embryo occurs already on a timescale td(j]/ e o) ~ 
to-\r P lr f ){a/R H )4Zt d . 

Planetesimals radially separated from the embryo by 
~ Rh pass the embryo's Hill sphere (and experience 
strong scattering) roughly every synodic period t syn cor- 
responding to the radial separation of Rh ■ 



t syn a n-\a/R H ). 



(37) 
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Requiring the damping time td to be shorter than t syn 
(which is necessary for planetesimals to stay in the shear- 
dominated regime) for embryos less massive than Mf is 
equivalent to demanding that 

V^{r p /r f )<a/R H or M e /Mf < {r p /rf)~ 3 . (38) 

Around more massive embryos with M e > Mf planetes- 
imal eccentricity is strongly damped between consecu- 
tive encounters with the embryo below its initial value if 

td(v/ e o) < tsyn Or if 



(r P /r f )(a/R H ) < a/R H 



or 



rp/r f < 1. (39) 



At the same time, eccentricity would not drop below 77 
prior to the next encounter if l|38l) is not fulfilled simul- 
taneously with (|39[1 . 

Comparing l|3*%|l and (|3*9"|) with the dispersion- 
dominated conditions (|3(J|) and l|32l) we see that depend- 
ing on M e /Mf and r p /rf there can only be two possible 
states in the system: (1) either planetesimals are scat- 
tered in a smooth, continuous fashion in the dispersion- 
dominated regime, with gas drag not capable of damping 
their eccentricities significantly between the consecutive 
approaches to the embryo, or (2) they are strongly scat- 
tered by the embryo in the shear-dominated regime at 
each approach and gas appreciably reduces their random 
velocities before the next encounter takes place. 

The discrete nature of planctesimal scattering in the 
shear-dominated regime is very important for determin- 
ing the approach velocity of planetesimals to the em- 
bryo. For example, if we were to assume that scatter- 
ing in the shear-dominated regime is continuous (like 
in the dispersion-dominated case) the average rate of 
eccentricity growth would have had the form de/dt ~ 
(Rh /a) /tsyn, since embryo increases planetesimal eccen- 
tricity by ~ Ruja every synodic period. Balancing this 
by the gas drag in the form 0, one would find the av- 
erage value of eccentricity to be ~ (Rh /a)(td/t syn ). It 
would however be a grave mistake to assume that this is 
the eccentricity with which planetesimal approaches the 
embryo. Indeed, it follows from l|35|l that the planetesi- 
mal eccentricity right before the encounter with the em- 
bryo is ~ (Rh/o) exp(— t syn /td) for M e < Mf, which is 
exponentially smaller than the average value of e! Thus, 
proper treatment of the shear-dominated regime taking 
the discrete nature of scattering into account is crucial 
for figuring out the initial conditions of the interaction 
process. This will have important ramifications for the 
question of accretion of these planetesimals as we demon- 
strate in 

3.2. Stokes drag. 

Planetesimals interact with the gas in the Stokes 
regime when r p > A and Re < Reb- The last condi- 
tion depends not only on the particle size r p but also on 
its velocity. We introduce another fiducial size, rg, which 
is defined as a planetesimal size for which Re = Reb at 
v = Av a : 
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Planetesimals with r p < rg always experience gas drag in 



the Stokes regime for M p < Mf (because post-scattering 
planetesimal velocity is < Av g , v. 
Re b ). 



Aw g , and Re < 



Embryos more massive than Mf endow shear- 
dominated planetesimals with velocity SIRh > Av g at 
each scattering episode. But even then planetesimals 
with sizes r p < rs satisfying condition 
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experience only the Stokes drag, i.e. their maximum ve- 
locity and physical size are never large enough for their 
Reynolds number to exceed Reb- Planetesimals scattered 
by more massive embryos would experience quadratic 
drag right after the encounter (even if only temporarily) 
before switching to the Stokes regime. From equation 
Ijl8(l one can easily see that damping between encounters 
is purely exponential for the Stokes drag independent of 
the planetesimal velocity (similar to the quadratic drag 
for e < 77, see H3 3.1(1 . 

Let us now turn to the dispersion-dominated regime. 
Balancing heating rate H25|) by the damping rate l|18fl 
with £ = (3/2)A/r p we find that 
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(42) 



Consequently, planetesimals are in the dispersion- 
dominated regime with respect to the embryo if 

M e/Mf > fe) (n»/*y)- 8 - («) 

Smaller planetesimals are in the shear-dominated regime 
and experience discrete scattering by the embryos. From 
equation Q42Jl we also find that v < Av g whenever 

M e /Mf < (^j ' {r p /r f )-\ (44) 

and that Re < Reb (and drag is in the Stokes regime) if 

M e /M f <(^j \r p /r f )-W (45) 

(see Figure |SJ|. The last equation is a dispersion- 
dominated analog of the condition (|41|l . 

3.3. Epstein drag. 

Smallest planetesimals with r p < A, are coupled to gas 
via the Epstein drag. Performing analysis analogous to 
that of H3 3. 21 one finds that in the dispersion-dominated 
regime scattering by the embryo maintains planetesimal 
random velocity at the level of 
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Planetesimals can only be in the dispersion-dominated 
regime with respect to the embryo if 

M e /M f > (^) 3 (r p /r f )- 3 (47) 
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Fig. 2. — Sketch of the planetesimal eccentricity (top) and in- 
clination (bottom) evolution due to the scattering by the embryo. 
Spikes correspond to close approaches to the embryo (which hap- 
pen at different intervals because of the changing semimajor axis 
separation). We assume that there is a weak continuous source 
of planetesimal excitation keeping random velocities above some 
minimum level (e m in an< i imin)- In reality behaviors of e and i are 
more erratic because of the scattering by distant embryos (see J1J. 



(see Figure • In this dynamical regime planetesimal 
velocity is below Av g only if 

M e /M f <(^y /2 (r p /r f )-^. (48) 

Planctesimals too small to satisfy (|47f) are in the shear- 
dominated regime and experience strong scattering by 
the embryo every synodic period, with their orbital ele- 
ments exponentially decaying between encounters (anal- 
ogous to the behavior in the case of Stokes drag, see 

3.4. Inclination evolution. 

It is easy to see from equations (II 61) and l |18| l that as 
long as gas drag is the only force affecting planetesimals 
after their encounter with embryo, planetesimal inclina- 
tion decays according to 

« We(*)/eo, (49) 

with eo,io being the post-encounter values of eccentric- 
ity and inclination and e(t) given by (|34[1 for quadratic 
gas drag. In the shear-dominated regime one might be 
tempted to think on the basis of (|49|l that i ^> e long 
after the encounter took place (when e <C eo). This ar- 
gument, however, assumes that io ~ eo and we now show 
that this is not the case. 

Indeed, even if io ~ Rh / a, inclination is exponentially 
small right before the next encounter (see an d |49j ) 
since the damping time is shorter than the synodic pe- 
riod in the shear-dominated regime. After this second 
encounter, eccentricity increases to ~ Rn/a, while the 



inclination does not go back to Rnja but, according to 
lf2I|. remains small, see Figure Subsequent action of 
gas drag prior to the next passage of the core further 
reduces planetesimal inclination, and so on. As a re- 
sult, we arrive at a very interesting conclusion: as long 
as planetesimals are in the shear-dominated regime their 
inclinations keeps decaying. Thus, all shear-dominated 
planetesimals "rain out" towards the disk midplane and 
collapse into geometrically thin layer. This emphasizes 
the importance of determining which part of the plan- 
etesimal population is shear-dominated with respect to 
the embryo. 

Thickness of this layer would be zero if embryo were 
on purely uninclined orbit and gas drag were the only 
process affecting planetesimal dynamics between encoun- 
ters. In reality, there are additional stirring agents which 
would keep the thickness of this subdisk finite (although 
still very small) and we consider them in 3H& 311 

3.5. Separation of different regimes. 

We now summarize what we have learned in H3 3.113 3.31 
about planetesimal dynamics in the proto-Solar nebula. 

First of all, it is clear from our previous discussion that 
the separation of different gas drag regimes sensitively 
depends on the relative values of fiducial planetesimal 
sizes A, r/, and rs given by the equations 10} , and 
<|40l) . In Figure|3]we display the scaling of these sizes with 
the distance from the Sun. The boundary between the 
quadratic and Stokes drag regimes is calculated assuming 
v w Av g and coincides with rs(a). We also display the 
curve r s to P {a) using equations (|13|I - H15(I — planetesimals 
with r p below this curve (shaded region) are tightly cou- 
pled to the gas and their dynamics cannot be described 
by equations ((TB|l - lfT^|) . 

From Figure |3 one can see that depending on the lo- 
cation in the nebula there could be three possible situa- 
tions: (a) when A < rs < rt, which holds for a < 3 AU 
in the MMSN, (b) when A < r $ < rs, which takes place 
for 3 AU < a < 9 AU, and (c) when rj < A < rs, which 
is valid for a > 9 AU. 

The first case pertains to the region of terrestrial plan- 
ets. We display possible regimes of planetesimal inter- 
action with gas in this part of the nebula in Figure 21 
They are classified according to the values of planetesi- 
mal size r p and embryo mass M e . Thick solid line sepa- 
rates planetesimals which interact with the embryo in the 
shear-dominated regime and experience discrete scatter- 
ing (hatched region to the left of the curve) from those 
which undergo the dispersion-dominated scattering (un- 
hatched region). In the terrestrial region, as we see from 
Figure this curve is set by equations (|30|l and (|32|l . 
Dashed curves denote the boundaries of different gas drag 
regimes: Epstein drag operates when r p < A; Stokes drag 
operates whenever Re < Re^ — in the shear-dominated 
regime this implies r p < rs for embryos that cannot kick 
planetesimals by more than Av g (M e < M f , see 
and restriction (|41|) for embryos that can 3 (M e > Mf). 
For even bigger planetesimals gas drag is quadratic and 

3 I n fa ct, even near embryos more massive than M e determined 
from 1411 planetesimals can spend some time in the Stokes regime: 
although they are in the quadratic drag regime right after the pas- 
sage of the embryo, their velocity rapidly decreases after that and 
their Reynolds number can drop below Rei, before the next en- 
counter takes place. 



7 




^ lil i i V"i' i iii t " r r" I 'Tt r~r it i i TO i"Vri hi 

10" 1 1 10 100 

Distance from the Sun a, AU 

Fig. 3. — Behavior of the length scales rf, rg, and A important 
for the planetesimal interaction with the gas (with different regimes 
marked on the plot), as a function of distance a from the Sun. 
Planetesimals with sizes in the shaded region have stopping time 
shorter than local Q~ 1 and are thus moving together with the gas. 
Complex shape of the boundary of this region is due to the different 
gas drag regimes. 
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Fig. 4. — Separation of planetesimals into the shear- and 
dispersion-dominated ones with respect to the embryo in the phase 
space of the embryo's mass M e and planetesimals size r p . This 
particular separation and numerical values indicated in the Figure 
pertain a = 1 AU (terrestrial planet region), where A < rg < r/. 
Thick solid line separates two dynamical regimes, with shear- 
dominated one being hatched. Thick dashed curves separate dif- 
ferent regimes of planetesimal interaction with the gas (explicitly 
indicated on the plot). See text for more details. 



planetesimals can be either in the shear- (small ones) 
or in the dispersion-dominated (large ones) regime. In 
the latter regime the dot-dashed line separates cases of 
equilibrium planetesimal random velocity being higher or 
lower than Av g , see equations (13 It and (|33|1 . When con- 
sidering the population of small planetesimals we should 
always keep in mind that our discussion in *13 3.113 372\ 
is valid only for planetesimals moving on almost Keple- 
rian orbits. Thus, a condition r p > r stop must be satisfied 
which restricts the validity of our results to planetesimals 
bigger than ~lm (see Figure EJ- Smaller planetesimals 
largely follow the motion of the gas. Embryo's mass is 
not allowed to exceed r]~ 3 / 2 Mf because more massive 
protoplanetary cores excite supersonic random velocities 
of surrounding planetesimals. 

As we have demonstrated in H3 3.41 small shear- 
dominated planetesimals sediment into the geometrically 
thin subdisk (thinner than the embryo's Hill radius) near 
the midplane of the nebula. We can readily see from 
the Figure 0] that in the terrestrial region this destiny 
awaits all planetesimals smaller than ss 100 — 200 m in 
size (less than ~ 10 13 — 10 14 g in mass) in the vicinity 
of embryos more massive than about 10 25 g. Near less 
massive embryos even larger planetesimals can belong 
to this dynamically cold population: only bodies bigger 
than about 1 km can escape this fate near 10 23 g mass 
embryo. Thus, gas drag can have important effect even 
on 0.1 — 1 km size planetesimals when the question of 
their interaction with the massive embryos is concerned. 

Figure0represents the separation of different gas drag 
and dynamical regimes at 5 AU from the Sun (A < Tf < 
rs), corresponding to the giant planet region (roughly 



the semimajor axis of Jupiter). Again, one should keep 
in mind that planetesimal smaller than about 1 m are 
tightly coupled to the gas at 5 AU, see Figure |3| Evi- 
dently, the structure of the M e — r p phase space is more 
complex in this part of the nebula. For example, plan- 
etesimals interacting with the embryo in the dispersion- 
dominated regime can now experience not only quadratic 
but also the Stokes drag. The reason for this is the lower 
Hill velocity £IRh for a given M e at 5 AU which re- 
duces the efficiency of planetesimal stirring by the em- 
bryo (see [23) and diminishes the equilibrium planetesi- 
mal velocity. As a result, Reynolds number can drop be- 
low Reb even for the dispersion-dominated planetesimals 
and drag can switch to the Stokes regime. Boundaries 
of different drag regimes are computed using equations 
l|4l"|) and (|43|) . Hatched region again represents M e and 
r p for which planetesimal scattering by the embryo oc- 
curs in the shear-dominated regime (its boundary is de- 
termined by equations [3Uj, |32|. and 03 )■ O ne can see 
that this happens for planetesimals that are somewhat 
smaller than at 1 AU: r p < 30 m for M e w 10 26 g and 
r p 100 m for M e ss 10 23 g. This is because gas density 
rapidly drops with the distance from the Sun diminish- 
ing the strength of dissipation, and this tendency cannot 
be counteracted by the longer synodic period at 5 AU 
(for a given M e ). Thus, somewhat less planetesimal ma- 
terial will be concentrated in the vertically thin subdisk 
of small bodies at 5 AU than at 1 AU, but the difference 
is not very pronounced. 

Finally, Figure displays the situation at 30 AU 
(roughly the semimajor axis of Neptune), in the region 
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Fig. 5. — Same as Figure Bl but for a = 5 AU (giant planet 
region), where A < rj < rg. 
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Fig. 6. — Same as Figure IH but for a = 30 AU (region of ice 
giants) , where r t < A < r g . 



of ice giants. Calculation of boundaries of different gas 
drag and dynamical regimes is performed using equa- 
tions P3|). gSJ-lgSJ, From Figure H we can easily 
see that the molecular mean free path A is very large in 
this part of the nebula making the Epstein gas drag im- 
portant for setting the boundary between the shear- and 
dispersion-dominated regimes. In this distant part of the 
nebula only bodies smaller than ~ 10 cm would be dy- 
namically coupled to the gas. Critical planetesimal size 
below which scattering is shear-dominated is 50 — 300 m 
for embryos with masses 10 23 — 10 27 g. This critical size 
is larger than at 5 AU because gas drag in the Epstein 
and Stokes regimes is more efficient than in the quadratic 
regime (for the same r p and M e ). 

3.6. Migration, gaps, and multiple embryos. 

Conservation of Jacobi constant ensures that any 
change in the epicyclic energy of planetesimal in the 
course of its scattering by the embryo is accompanied by 
the change in planetesimal semimajor axis. As a result, 
planetesimal surface density distribution gets perturbed 
by the embryo and a gap might form (Ida & Makino 
1993). For the shear-dominated planetesimals gap open- 
ing can be especially important, because in this dynami- 
cal regime planetesimal guiding centers are moved away 
from the embryo's orbit by ~ Rh in a single passage at 
h w Rh, see (|2*3")) (Rafikov 2001). If this happens, ac- 
cretion of planetesimals by the embryo can be severely 
affected (Rafikov 2003a). 

At the same time, gas drag causes orbital decay of 
planetesimals between encounters — they migrate to- 
wards the Sun. This drift moves planetesimals which 
are located inside of the embryo's orbit further from it, 
facilitating the gap opening. On the other hand, on the 
outer side of the embryo's orbit gas drag causes plan- 
etesimals to drift towards the embryo, and this tends to 



oppose the gap formation. As long as planetesimals are 
less massive than the embryo, they are all repelled by 
scattering in the same way independent of their masses. 
Thus, the question of whether the gap outside of the em- 
bryo's orbit is cleared or not depends only on the damp- 
ing timescale; there exists a critical size r m i g such that 
planetesimals of this size initially at h ~ Rh from the 
embryo, after being repelled by roughly Rh, can migrate 
back the same distance in a synodic period. Only these 
planetesimals would be accreted by an isolated embryo: 
planetesimals with r p > r mig are too big to be brought 
back to the embryo by the gas drag, and a gap forms 
preventing their further accretion. Planetesimals smaller 
than r m i g migrate so fast that in a synodic period they 
cross the embryo's orbit and are lost to the inner disk. 
Embryo's accretion would then be rather inefficient as 
well. 

These problems arise only for an isolated protoplane- 
tary core. However, during the intermediate stages of 
planet formation as a consequence of oligarchic growth 
(Kokubo & Ida 1998) there would be many embryos 
present in a disk at the same time. Their orbits should 
not be very widely separated: even if it were the case ini- 
tially, subsequent increase in the embryo masses caused 
by the accretion of planetesimals would make these sep- 
arations not too large (see below) in terms of their Hill 
radii (because Rh expands as embryo's mass increases, 
see |20| ) and this is what is important for the dynamics. 

When such a "crowded" population of protoplanetary 
cores (which is a natural outcome of oligarchic growth) 
is present in the disk, gap formation is no longer an is- 
sue: although a particular embryo repels planetesimals 
and tends to open a gap, scattering by another nearby 
embryo pushes planetesimals back and spatially homog- 
enizes them before they approach the first embryo again. 
Gap formation is thus suppressed and accretion can pro- 
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ceed almost uninhibited. On the other side of the prob- 
lem, although small shear-dominated planetesimals mi- 
grate through the orbit of some particular embryo be- 
cause of the gas drag, there are many other embryos 
in the inner disk which scatter planetesimals back and 
forth, so that their inward drift looks more like a ran- 
dom walk through the nebula. In the course of such 
a diffusion through the protoplanetary disk planetesi- 
mals have a high chance of being accreted by one of the 
many embryos (as long as planetesimals stay in the shear- 
dominated regime, see IJB}. 

At a glance, it seems improbable that embryos can 
remain on purely circular and uninclined orbits in a 
"crowded" configuration because at small radial separa- 
tions they would strongly scatter each other and very 
quickly excite large random velocities. However, one 
should remember that in the course of oligarchic growth 
most of the solid mass is locked up in planetesimals and 
not in embryos. This is true until M e reaches the isola- 
tion mass Mi SO defined as 

M lso « Mq (4nZ p a 2 /M f /2 « 6 x 10 26 g a 3 /* (50) 

if embryos grow predominantly by accretion of shear- 
dominated planetesimals and 



M, 



6 x 10 
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) (47r£ p a 2 /M ( 
g (V/ilR H f 2 a 3 /* 



3 , 3/2 (v/sir h . 



3/2 



(51) 

if embryo growth is dominated by accretion of large 
dispersion-dominated planetesimals with typical random 
velocity V (e.g. see Rafikov 2003b). Massive planetes- 
imal population exerts dynamical friction on embryos 
transferring random energy of their epicyclic motion to 
planetesimals, which keeps embryo eccentricities small 
(see Kokubo & Ida 1995). Planetesimal velocities, in 
turn, are damped by the gas drag 4 which allows dynam- 
ical friction to continue being effective. 

Planetesimal disks have presumably contained both 
massive dispersion-dominated planetesimals and small 
shear-dominated bodies at the same time. We will as- 
sume that the latter comprise some fraction % < 1 of 
the total planetesimal surface density S p . One can show 
(Stewart k, Ida 2000; Rafikov 2003c) that these two pop- 
ulations of planetesimals damp random velocity v e of pro- 
toplanetary cores at a rate 
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where second term represents the contribution of large 
dispersion-dominated planetesimals (we are dropping 
here all constant factors and the Coulomb logarithm), 
while the first term is due to the small shear-dominated 
bodies 5 . Evidently, small bodies are more important 

4 In the end, the energy of the embryo's epicyclic motion gets 
damped by the gas drag, with planetesimal population acting as 
an intermediary. 

5 Effect of dynamical friction by the shear-dominated bodies can 
be obtained from that caused by the dispersion-dominated plan- 
etesimals by setting v = VIRh and reducing planetesimal surface 
density by a factor of \. 



for "cooling" the cores than large planetesimals provided 
that 



X 



> 



QR 



H 



V 



(53) 



Note that because V > HRh this inequality can be 
fulfilled for rather small x ( e -g- f° r X = 0.02 for 
V = 3QRh)- In the case opposite to (|53() dynamical 
friction on the protoplanetary cores is controlled by the 
large dispersion-dominated planetesimals. 

We demonstrate later in that the relative impor- 
tance of shear- or dispersion-dominated populations for 
the growth of mass of protoplanetary cores is determined 
by a condition different from (|53|l . Namely, protoplane- 
tary cores grow predominantly through the accretion of 
small shear-dominated bodies whenever 



x> P 1/2 m H /vy 



(54) 



where parameter p <C 1 is defined by equation (|78fl . In 
this accretion regime one can also demonstrate that ra- 
dial separation of embryos h m Rh and their gravita- 
tional scattering has a discrete character: random veloc- 
ities of embryos are strongly excited during their close 
approaches to each other, but dynamical friction tends 
to damp them before the next encounter occurs; this is 
similar to scattering of the shear-dominated planetesi- 
mals in the presence of gas drag (see 3.113 3~1|) . A 
particular embryo grows mainly by accretion of material 
(small planetesimals) from the annulus around its orbit 
(the so called "feeding zone" ) with the radial width « Rh 
leading to isolation mass in the form (|50|l . For a system 
of embryos to remain most of the time on uninclined 6 
and roughly circular orbits with radial spacing h « Rh, 
dynamical friction timescale tdf must be shorter than the 
average time between the embryo encounters t syn given 
by (|37|l . Estimating tdf from (|52|l one can find that this 
is possible only for embryos with masses satisfying 



x + (i-x) 
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(55) 



Note that if the whole planetesimal disk is shear- 
dominated, i.e. x = 1) then embryos can stably stay on 
closely packed (h ss Rh) almost circular and uninclined 
orbits all the way until they reach the isolation mass 
(modulo constant factors), compare with (|50J) . However, 
if x < 1 critical embryo mass goes down compared to 
M 

Whenever the amount of solid mass in shear-dominated 
planetesimals is so small that (|54f) is violated and % ^ 
p x / 2 (Oiiij /V) 2 , core growth is determined mainly by the 
accretion of large dispersion-dominated planetesimals. 
One can demonstrate in this regime that as a result 
of oligarchic growth radial separations of protoplanetary 
cores become » V/Q (Ida & Makino 1993; Kokubo & 
Ida 1998) which is larger than Rh- Feeding zones of 
embryos are larger than in the shear-dominated case: 

6 Analogous to the planetesimal case, the shear-dominated scat- 
tering of embryos by embryos is very inefficient in exciting vertical 
velocities. 
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dispersion-dominated planetesimals can be accreted from 
within V/ft in the radial direction (i.e. the width of the 
feeding zone is again roughly equal to the radial separa- 
tion between the embryo orbits). Since orbits of embryos 
are well separated (h V/ft > Rh) their gravitational 
scattering is not so dramatic as in the case correspond- 
ing to H54fl and eccentricity growth has a character of 
a random walk. Using equation lll'L'1) and the fact that 
embryos encounter their neighbors roughly once every 
tsyr^Rn/h) w t syn (Q,R H /V) (where t syn is defined by 



) we find growth rate of random velocity to be 



dt 
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(56) 



For the system of embryos to be dynamically stable, the 
equilibrium value of random velocity v e (obtained by bal- 
ancing scattering rate |56| with the damping due to the 
dynamical friction [52]) has to be less than V — orbits 
of cores should not cross. One can easily find that this 
is possible only provided that 

3 / 2 / T/ \ 3 / 2 
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(57) 



If planetesimal disk contains only dispersion-dominated 
planetesimals, i.e. x — 0, then population of embryos 
can be dynamical stable only until embryos reach cor- 
responding Mi SO given by (|51|) . In Figure we display 
different regimes of dynamical friction and planetesimal 
accretion as a function of \ and V/ClRff. 

Thus, whenever planetesimal disk is one-component 
(i.e. has either only shear-dominated or only dispersion- 
dominated planetesimals), it is the isolation mass that 
sets a limit on the maximum core mass below which pro- 
toplanetary cores co-exist on closely spaced but dynam- 
ically cold orbits. However in the more realistic case of 
a disk consisting of both populations (i.e. < \ < 1) 
dynamical stability is violated at masses smaller (some- 
times much smaller) than Mi SOl as equations l|55|) and 
l|57l) clearly demonstrate. For instance, let us con- 
sider the region of giant planets at a = 5 AU where 
p 10~ 3 according to l|78|) . Assuming that most of the 
surface density is in large, 1-km size planetesimals, we 
find their random velocity V ~ 2.2flRn to be set by 
quadratic gas drag Then, whenever the mass frac- 

tion of small shear-dominated planetesimals forming al- 
most two-dimensional subdisk is \ S 6 x 10 -3 , both core 
accretion rate and dynamical friction are set only by large 
planetesimals. For such small x t ne system of embryos 
is dynamically stable until M e reaches the isolation mass 
JHIJl, see (jnZJ). For 6 x 10~ 3 < \ < 4 x 10 -2 situation is 
different: dynamical friction acting on cores is still deter- 
mined by large planetesimals but core accretion rate is 
now set by the abundance of small bodies. Orbits of em- 
bryos accreting small bodies tend to get closely packed 7 

7 When x passes through the value p 1 / 2 (£IRh /V) 2 equilib- 
rium radial separation between the embryo orbits discontinuously 
changes from V/Vl to Rh- 
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Fig. 7. — Different regimes of core accretion and dynamical fric- 
tion in the presence of two populations of planetesimals — shear- 
and dispersion-dominated. Separation of regimes is presented in 
terms of V/QRh — velocity of large planetesimals scaled by the 
embryo Hill velocity and \ — planetesimal mass fraction locked up 
in s mall bodi es. Solid and dashed lines correspond to conditions 
1531 and I5H . 



and this increases gravitational scattering between them. 
As a result, dynamical instability sets in at a consid- 
erably smaller mass M ecr « Mi SO (QRH /V) 6 , where 
Mi SO is now given by l|50l) : for V = 2.2VLRh one finds 
M e . cr sa 0.01M lso . Finally, when x > 4 x 10~ 2 both dy- 
namical friction and core accretion rates are determined 
by small shear-dominated bodies and M e cr w M iso x 3 ^ 2 , 
(M iso is again given by jHOl); oni y when x = 1 can insta- 
bility be postponed until M e reaches Mi SO , 

From this example one can see that the situation in 
two-component planetesimal disks is very different from 
that in one-component disks. Even small admixture 
of shear-dominated planetesimals can completely change 
the dynamics of a system of protoplanetary cores. Thus, 
it becomes especially important to know the exact value 
of x and to follow its evolution in time. Table 1 briefly 
summarizes our findings by delineating the conditions 
under which shear-dominated or dispersion-dominated 
planetesimals control core accretion and dynamical fric- 
tion. 

As cores grow, their Hill radii increase and feeding 
zones overlap leading to occasional mergers of cores. This 
keeps their orbital separations from becoming too small 
but does not lead to dynamical instability because of 
the dynamical friction. Thus, the overall picture of oli- 
garchic growth described above is not affected by mergers 
of embryos. However, as soon as M e increases beyond 
the threshold given by (|55(l or (|57(l . dynamical friction 
can no longer keep embryos on kincmatically cold orbits, 
their eccentricities and inclinations start to grow, and 
cores finally switch into the dispersion-dominated regime 
with respect to each other. Similar effect has been ob- 
served by Kokubo & Ida (1995) in N-body simulations, 



11 



although they were dealing with the gas-free environment 
which made dynamical friction less effective since plan- 
etesimals were dynamically hot. In this work we are not 
going to follow this more violent stage of planetesimal 
disk evolution. 

4. LOWER LIMIT ON PLANETESIMAL VELOCITY. 

We have demonstrated in 3.1l that planetesimal ran- 
dom velocity in the shear-dominated case rapidly decays 
as a result of gas drag after scattering by the embryo. If 
the damping time td is much shorter than the time be- 
tween encounters, planetesimal velocity right before the 
next approach to the embryo would essentially be zero 
[oc exp(— tgyn/td) to be more exact]. As a result, any 
subdominant sources of stirring which would normally be 
negligible become important in maintaining random mo- 
tions at a finite value. These effects then determine the 
floor below which planetesimal random velocity cannot 
drop. Here we identify two such effects — scattering by 
the distant embryos and stirring by the large dispersion- 
dominated planetesimals — which determine the mini- 
mum horizontal and vertical random velocities of plan- 
etesimals respectively. We consider these two processes 
separately. 

As we demonstrated in i j3 3.61 proto-Solar nebula 
should contain a population of dynamically cold em- 
bryos separated by roughly Rh/N — we parametrize 
the uncertainty in the radial separation of cores by the 
number of embryos N per Rh in radius; if cores grow 
mainly through accretion of shear-dominated planetes- 
imals N « 1, while if they mostly accrete dispersion- 
dominated bodies N « (flR H /V) < 1, see %TTH1 Sur- 
face number density of embryos is w N/(2naRH)- 

Distant embryos scatter a given planetesimal quite 
frequently because both the velocity of incoming em- 
bryos (which is determined by the shear in the disk) and 
their number increase linearly with the radial separation 
h. As a result, although scattering by the nearest em- 
bryos is discrete, scattering by the cores more distant 
than some critical h c should be considered as a contin- 
uous process, similar to the scattering in the dispersion- 
dominated regime. To determine h c we notice that scat- 
tering switches from the discrete to continuous mode 
when the average time between the passages of embryos 
separated from a given planetesimal by less than h c be- 
comes shorter than the typical timescale on which plan- 
etesimal velocity would evolve otherwise. For planetes- 
imals with r p > r st0 p this typical timescale is the gas 
damping timescale 8 td- 

The rate T(h c ) at which embryos with \h\ < h c pass a 
particular planetesimal due to the shear in the disk is 



is given by T{h c )td ~ 1, meaning that 
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Boundary between the discrete and continuous scattering 

8 Note that for planetesimals smaller than r s top gas damping 
time is shorter than Q~ 1 . However the duration of the gravitational 
interaction with the distant embryo is always ~ f2 ~ 1 meaning that 
the typical velocity evolution timescale for planetesimals with r p < 
"Tstop is not t,j but rather Cl~ 1 . This, however, is only important 
for very small planetesimals, r p < 1 m (see Figure Uj, which are 
not covered by this study anyway. 
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Thus, h c > Rh for TV < 1 and t d < t syn . 

Scattering by the embryos with \h\ 3> h c occurs so 
frequently compared to the gas damping timescale that 
the random component of scattering averages to zero and 
eccentricity stirring is given by (|22() . We find that 
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Balancing this stirring by the gas damping de 2 /dt — 
—e 2 /td and using (|59() we find the equilibrium value of 
eccentricity 



3N Rh td Rh 



4-7T a t* 



Rh 
h, 



~~r~ ) ■ (61) 



Using we evaluate 
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for quadratic, Stokes, and Epstein drag regimes respec- 
tively. These estimates imply that the minimum horizon- 
tal random velocities of small planetesimals are below the 
Hill velocity, as they should be. For 10 m planetesimals 
stirred by 10 25 g embryos e m in corresponds to velocities 
of the order of 1ms" 1 in the inner part of the proto-Solar 
nebula, dropping to w 0.5 m s -1 at 30 AU. 

At the same time, scattering by distant embryos can- 
not maintain the inclinations of small planetesimals at 
a finite level. Excitation of the vertical velocity by an 
encounter with an embryo separated even by Rh from 
planetesimal is weakened compared to the excitation of 
horizontal velocity by the geometric factor ia/Rn <C 1. 
As a result, the growth rate of inclination due to the em- 
bryo scattering scales as i 2 (exactly like gas drag) and for 
td ^$ tsy n S as drag unconditionally dominates. We now 
consider if stirring by planetesimals (and not embryos) 
can keep nonzero inclinations of small bodies. 

Gas drag acting on planetesimals bigger than 0.1-1 km 
(depending on the location in the nebula, see WS 3.5f) is 
too weak to prevent them from staying in the dispersion- 
dominated regime with respect to embryos. Gravita- 
tional interaction of these massive planetesimals with 
small bodies is certain to take place in the dynamically 
"hot" regime (because for the same physical velocity the 
Hill radius for the planetesimal-planetesimal scattering 
is much smaller than Rh for the embryo-planetesimal 
scattering) . 

Rafikov (2003c) has demonstrated that velocity ex- 
citation by planetesimals sensitively depends on the 
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Fig. 8. — Schematic representation of the mass spectrum of solid 
bodies in the disk (thick solid curve). Three distinct parts of the 
spectrum can be singled out: small bodies (dAf /dm oc m~@ , < 
j3 < 2, most of the mass is in biggest objects), merger spectrum 
(dAf I dm oc m~ a , 2 < a < 3, most of the mass is in smallest 
objects), and a set of embryos detached from continuous mass dis- 
tribution (see J1J. Arrows describe different mass transfer routes 
between the components of the mass spectrum: solid arrow is for 
the direct accretion of large (dispersion-dominated) planetesimals 
by embryos, while dashed arrows describe fragmentation of plan- 
etesimals with further shear-dominated accretion by embryos. See 
J7|for more details. 



planetesimal mass spectrum. For a given differen- 
tial surface number density distribution of planetesimal 
masses dAf /dm the inclination stirring can be written as 
(Rafikov 2003c) 



dt \ V 



ria\ 2 a 2 
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dm 



(65) 



where V is the random velocity of large dispersion- 
dominated planetesimals (which wc for simplicity set 
constant, independent of m). Shear-dominated planetes- 
imals are not efficient at velocity excitation (but they 
matter for dynamical friction). Because of the stirring 
by embryos V should be some multiple of QRh , but ow- 
ing to the action of gas drag (and planetesimal dynamical 
friction) it is not higher than several flRn, see equations 
123, 02>, and (2J- 

Numerical simulations very often produce planetesimal 
mass spectra such that dAf /dm oc m~ a with a sa 2.5 
within a wide range of masses (Kokubo & Ida 1996). In 
a disk with such spectrum most of the mass is concen- 
trated at the lower end of distribution while most of the 
stirring is done by its upper end (Rafikov 2003c) . We will 
assume that planetesimal mass spectrum has this form 
for mo < m p < m*, where we somewhat arbitrarily take 
mo = 10 16 g and m* sa 10 22 g (roughly 1-km and 100-km 
size planetesimals); see Figure |H1 for a schematic picture 
of the assumed size distribution. Excitation of inclina- 
tion by the largest planetesimals (m p w to*, which are 



still much smaller than the embryos) can be expressed in 
terms of m , m*, and planetesimal surface mass density 
S p (dominated by planetesimals with m p « mo) roughly 
as (Rafikov 2003c) 
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In writing down this expression we have assumed that 
large planetesimals are numerous enough to give rise to 
a continuous rather than discrete stirring. Balancing l|66|l 
with —i 2 /td we find the equilibrium value of inclination: 
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Evaluating this expression for V = 3HRh, M e = 10 
E p given by (J5J, tj given by lfTU|l - lfT2l . and our adopted 
values of mo, m*, and a = 2.5 we find that 
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(68) 
(69) 
(70) 



for corresponding gas drag regimes. These values of in- 
clination correspond to very small vertical velocities of 
planetesimals, ~ 1 cm s _1 in the vicinity of 10 25 g proto- 
planetary cores at 30 AU. Evidently, since i m in <C e rn i n 
(see jn2]-jnil) this type of stirring would not be able to 
affect minimum eccentricities of planetesimals. 

Scattering by large planetesimals would vertically per- 
turb not only small bodies but also embryos. Balancing 
the stirring rate 1)66(1 with the dynamical friction rate 
Q52[l. we find that the minimum inclination of the em- 
bryo i e ,min IS about 



Rh ( V 



\nR 



H 



m 

Me 



m 



3-Q 



1/2 



V 



QR 



H 



10 



-3 r h 



V 



-1/2 



( M e 



nR H J V.10 25 g 



-1/2 



i + x 



( V 



\flR H 



provided that the time between encounters with large 
planetesimals (which dominate vertical stirring) is 
shorter than tdf (for continuous approximation to hold) . 
As a bottom line, we may conclude that scattering 
by large planetesimals keeps inclinations of small bod- 
ies relative to protoplanetary cores at the level of ~ 
10~ 3 (Rh /a): most likely via the stirring of cores. 

There are other possibilities for maintaining inclina- 
tions of shear-dominated planetesimals at some mini- 
mum level. One of them is the gravitational scattering 
between small bodies themselves which transfers the en- 
ergy of random motion from horizontal into the vertical 
direction. This is likely not to be important because dy- 
namical relaxation of small bodies is very slow. Another 
possibility is a gravitational instability in the thin layer 
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(P. Goldreich, private communication) which may excite 
random velocities of small constituent bodies. Assuming 
that a fraction \ of solid mass is locked up in small shear- 
dominated bodies, one can find that instability would 
operate if the planetesimal velocity dispersion is below 
X7rGS p /J7 ~ 30 cm s _1 x (usual Toomre criterion), and 
this is still below e m i n VLa even for x = 1. Finally, in- 
elastic collisions between small planetesimals can heat 
the disk vertically, and we elaborate more on this in the 
next section. 

5. INELASTIC COLLISIONS. 

Inelastic collisions between planetesimals act as effi- 
cient source of damping. The escape speed from the sur- 
face of a 100 m body is about 10 cm s _1 and planetesi- 
mals are typically moving with higher velocities (see |G2j - 
|64j ) . which means that (a) they would lose a lot of energy 
in high-energy collisions, and (b) gravitational focusing is 
unimportant and collision cross-section is almost equal to 
the geometric cross-section of colliding bodies. Assum- 
ing that planetesimals lose a fraction ~ 1 of their energy 
when they collide, we estimate the velocity damping in 
inelastic collisions to be given by 



dt 



2 X S P V h 
PpTp V z ' 



(72) 



where we have again assumed that only a fraction x of 
solid mass is in the shear-dominated planetesimals. From 
this expression it is clear that inelastic collisions lead to 
exponential damping of velocity (if Vh/v z ~ const) on a 
timcscale 
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Importance of inelastic collisions is judged by compar- 
ing U ne i with td — d amp ing time due to the gas drag. 
Using equations (|10fl - l|12fl we find 



td/U n ei ~ 0.25 x (vh/vz) 

td/Unel ~ 0.1 X( v h/V z ) 
td/Unel » 0.01 X (Vh/Vz) , 



(74) 
(75) 
(76) 

for quadratic, Stokes, and Epstein gas drag regimes. 
These estimates clearly demonstrate that td ;$ Uriel 
whenever v z ~ Vh (i ~ e), meaning that inelastic col- 
lisions are unimportant when planetesimal velocities are 
roughly isotropic. This is always true in the dispersion- 
dominated regime allowing us to neglect inelastic colli- 
sions in this case. 

Situation is different in the shear-dominated regime 
where one can easily have v z <C Vh (see HJ. Let's assume 
that planetesimal inclination right before the encounter 
with some embryo is very low. Right after the scatter- 
ing event inclination remains roughly the same (see |24j ) 
while eccentricity goes up to ~ (Rjj/a) 3> i m in] as a con- 
sequence, ti ne i immediately becomes much shorter than 
td and inelastic collisions suddenly become more impor- 
tant than gas drag. However, this is a very transient 
stage because inelastic collisions not only dissipate en- 
ergy but also isotropize the velocities of colliding bodies. 
If planetesimals were colliding like rigid balls, bouncing 



off after collision, one would expect roughly isotropic re- 
coil velocities. Consequently, one would expect v z rs vh 
after every planetesimal has experienced a single physical 
collision, i.e. after time i,- ne i has passed since the scatter- 
ing by the embryo. This immediately reduces the impor- 
tance of inelastic collisions and makes gas drag more im- 
portant again for the random velocity dissipation shortly 
after scattering has occurred, see lJTI |l -lf77) )l . 

After that, according to 149(1 . inclination decays slower 
than eccentricity does meaning that the role of inelas- 
tic collisions keeps decreasing. If the time interval be- 
tween the successive encounters with the embryo is long 
enough, the eccentricity decay will stop at the asymp- 
totic value e m i n (of course, with occasional oscillations 
due to scattering by distant embryos) while the incli- 
nation would continue to decay further until it reaches 
imin "C e m i n . At this stage v z /vh goes down and physical 
collisions can again start occurring quite frequently; how- 
ever, the imminent isotropization of velocities after ev- 
ery collision limits their importance (in comparison with 
that of the gas drag) only to short periods of time. As 
a result, the timescale of velocity damping between en- 
counters with embryos should still be very close to the 
gas damping time scale td- 

Physical collisions might affect the determination of 
imin since they can be effective at pumping the energy 
of horizontal motions into vertical ones. The minimum 
velocity anisotropy v z /vh they provide can be determined 
from the condition td/ti ne i ~ 1 and turns out to be 
roughly 0.25%, 0.1%, and 0.01% for typical parameters of 
quadratic, Stokes, and Epstein drag regimes, see equa- 
tions (|74|I - H76|I . Scattering by distant embryos maintains 
horizontal random velocities of small planetesimals at the 
level given by l|fi^|l - (|Mjl . and this degree of anisotropy re- 
sults in i m in ~ (v z /vh) e min which for 10 m planetesimals 
is about 0.01x(Rh/cl) at 1 AU and 5 x 10~ 4 x(iW a ) at 
30 AU. Comparing this with fi fy - 1(70 )1 and JE) we may 
conclude that for not too small values of x 0.1) in- 
elastic collisions are important for setting the minimum 
value of relative embryo-planetesimal inclination in the 
presence of gas drag in quadratic or Stokes regimes (in 
the inner parts of the proto-Solar nebula, inside 5 — 10 
AU). Whenever shear-dominated planetesimals are af- 
fected by the gas drag in the Epstein regime, inelastic 
collisions can compete with stirring by large planetesi- 
mals only for x ~ 1- 

This discussion has assumed rather idealized model 
of planetesimal collisions (rigid balls). In reality high- 
energy impacts are likely to be catastrophic, leading to 
the disruption of participant bodies. Then the amount of 
kinetic energy transferred into vertical motions and the 
degree of isotropization would be determined by the ejec- 
tion velocities of resulting debris. Observations of veloc- 
ity dispersions in collisional families of asteroids (Zappala 
et al. 1996) and results of numerical simulations (Michel 
et al. 2003) indicate that ejection velocities are consider- 
ably smaller than the initial planetesimal velocities. This 
slows down the velocity isotropization compared to the 
case of collisions of "hard balls" . Nevertheless, we do not 
expect this to seriously change the general picture out- 
lined before. Besides, it is not at all clear how anisotropic 
would ejection velocities be in the case of planetesimals 
(10-100 m in size, strength-dominated fragmentation) 
colliding at several tens of m s _1 , since collisions between 
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10-km asteroids (gravity-dominated fragmentation), on 
which information exists, occur at relative speeds of sev- 
eral km s" 1 . 

6. ACCRETION OF LOW-ENERGY PLANETESIMALS. 

Protoplanetary embryos grow by accreting planetesi- 
mals. Accretion in differentially rotating disks is intrin- 
sically complicated because of the three-body gravita- 
tional interaction of the two merging bodies and the cen- 
tral mass. One can however greatly simplify this prob- 
lem by treating accretion as a two-body process while 
approximately taking into account three-body effects by 
limiting from below planctesimal approach velocity rela- 
tive to the embryo by QRh ■ This is a direct consequence 
of the shear present in the disk. 

With this in mind we can write a rather general for- 
mula for the protoplanetary mass accretion rate: 
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n p m p R 2 e v a ( 1 + 
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(77) 



Jp H Rhv z \ v a 
where R e is the physical size of the embryo and the last 
equality holds for v a <C v esc . Dimensionless physical size 
of the core p relative to its Hill radius is independent of 
the core mass but varies with the distance from the Sun: 
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Clearly, R e <g; Rh and p <C 1. 

Scattering of planetesimals by the embryo in the 
dispersion-dominated regime always tends to maintain 
v z ps v ps v a > £IRh, thus 
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dt \ v 

In the shear-dominated regime approach velocity v a is 
almost independent of planetesimal eccentricity or incli- 
nation and is ps QRh- At the same time, vertical velocity 
of planetesimals v z setting the thickness of planetesimal 
disk is smaller than QRh, and we find that 
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dt ' ' v z 

Derivation of equation l|77|) has implicitly assumed 
that the maximum impact parameter at infinity with 
which planetesimal can be accreted by the embryo 
R e \/1 + Ve SC /v% is smaller than the planetesimal disk 
thickness v z /ft. In the shear-dominated case, when 
v a ~ QRh, this sets a limit on the vertical velocity of 
planetesimals at which their accretion can still be de- 
scribed by equation (f77|l : 

v z > v z , cr = P 1/2 QR H ps 0.07 QR H a^j 2 (81) 

(Greenberg et al. 1991; Dones & Tremaine 1993). Equa- 
tion (I80|) is applicable only for v z > v z _ cr . Whenever 
v z < w z , cr planetesimal disk is very thin and embryo can 
accrete the whole vertical column of material it encoun- 
ters (eccentricity is only restricted to be smaller than 
Rh /a) ■ In this case one can easily show that 
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(82) 



This is the highest possible accretion rate of planetesi- 
mals by the protoplanetary core that can be achieved in 
the planetesimal disk. 

It is possible that most of the planetesimal mass has 
been concentrated in bodies with sizes of about 1 — 10 
km (merger products in Figure |HJ which interact with 
embryos in the dispersion-dominated regime. According 
to l|79fl the growth timescale of the embryo due to the 
accretion of these bodies is 
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Using equations (|27Jl , , and (|4^|l we find that the time 
needed to build 10 25 g protoplanetary core by accretion 
of such dispersion-dominated planetesimals is ~ 10 5 yr 
at 1 AU, w 4 x 10 6 yr at 5 AU, and ps 3 x 10 7 yr at 30 
AU. 

At the same time, as we hypothesized in 3.61 a frac- 
tion x < 1 of planetesimal mass could have been in small 
bodies which are shear-dominated with respect to em- 
bryos. This mass is concentrated in a vertically thin 
population and comparing l|55 )l -l(7D |l & with (J&TJ we 
find vertical velocity of planetesimals in the subdisk to 
be smaller than v z ^ cr . Thus, small shear-dominated bod- 
ies can be very efficiently consumed by embryos, and we 
find using l|H2^l that the embryo's growth timescale due 
to their accretion is 
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This is w lO 3 ^ 1 yr at 1 AU, w 10 4 x _1 yr at 5 AU, and 
« 2 x 10 5 x _1 yr at 30 AU. Thus, even if only about 1% 
of the mass in solids is locked up in the shear-dominated 
planetesimals (x = 10 -2 ) they would still dominate the 
accretion rate of the embryo because tg tS d < t g ^- Com- 
paring 179|) and (|82|l with E p lowered by a factor of \ 
we find that accretion of small shear-dominated bodies 
from very thin disk is more important for embryo growth 
than accretion of large dispersion-dominated planetesi- 
mals whenever x is such that condition (|54|l is fulfilled. 

7. DISCUSSION. 

The outlined picture of planetesimal dynamics in the 
gaseous nebula naturally divides planetesimals of differ- 
ent sizes into two well-defined populations with respect to 
gravitational scattering by a set of protoplanetary cores. 
One is a "hot" population of bodies bigger than 0.1—1 km 
which interact with embryos in the dispersion-dominated 
regime and have large inclinations so that their scattering 
is an intrinsically three-dimensional process. The other 
is a "cold" population of smaller planetesimals (sizes be- 
low 0.1-1 km) which interact with embryos in the shear- 
dominated regime. These planetesimals tend to be con- 
fined by the action of the gas drag to a vertically thin 
disk [with a thickness of ~ 10~ 3 Rh] near the nebular 
midplane. 

Difference between these populations is most striking 
when the accretion of planetesimals by the embryos is 
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concerned. Accretion of dynamically hot planetesimals 
allows 10 25 g embryo to double its mass in ~ 10 5 yr in 
the terrestrial planet region and in several 10 7 yr in the 
region of ice giants. The latter timescale is quite long 
from the cosmogonical point of view. At the same time, 
accretion of cold planetesimals is about 100 times faster, 
with mass doubling timescale of ~ 10 3 yr and ~ 10 5 yr in 
the inner and outer parts of the nebula, provided that all 
solid mass is locked in these small planetesimals. Forma- 
tion of gaseous atmospheres around massive protoplane- 
tary cores further accelerates accretion of small bodies: 
gaseous envelopes can be very efficient at trapping small 
planetesimals (Inaba & Ikoma 2003) , thus increasing the 
capture radius R e and shortening the growth timescale. 
Protoplanetary cores would also grow by merging with 
other cores (since in our picture they reside on closely 
spaced orbits as long as M e < M e ^ rit , see M3 3.6[) : appar- 
ently, the faster cores grow, the closer their orbits are in 
the Hill coordinates, and the faster they merge. 

It is possible that most of the solid mass in the proto- 
Solar nebula was concentrated from the very start in 
small (10— 100 m) dynamically cold bodies. In this case 
X ~ 1 and growth of embryos should be very fast, see 
equation It is, however, equally possible that most 
of the mass was initially locked up in large (1 — 10 km) 
dispersion-dominated planetesimals and not in the dy- 
namically cold population of small bodies, i.e. \ !• 
In this case, although there is a huge reservoir of solid 
mass potentially available for the accretion, embryos can 
hardly make use of it because hot planetesimals are ac- 
creted rather inefficiently, see equation l|83l) . At the same 
time, cold population which can potentially allow a vig- 
orous growth of the embryo might simply not contain 
enough surface density to ensure high enough accretion 
rate (accretion timescale is inversely proportional to the 
fraction of mass contained in small planetesimals, see 
|84j ) . Thus, unless enough mass (x > 10 -2 ) is trans- 
ferred from hot to cold planetesimals, embryos would 
grow slowly accreting large dispersion-dominated plan- 
etesimals. 

A natural process for transferring mass from the big 
bodies into small ones is planetesimal fragmentation (see 
Figure |HJ), which should naturally be taking place in the 
presence of massive embryos. Indeed, the escape speed 
from the surface of 1 km body is about 1 m s _1 , while 
the Hill velocity of 10 25 g embryo is w 50 m s _1 at 1 
AU and is « 10 m s _1 at 30 AU. Thus, 1-km plan- 
etesimals possibly containing most of the solid mass in 
the disk would collide with kinetic energy far exceeding 
their gravitational binding energy. Depending on their 
internal strength, parent bodies can be disrupted into a 
number of smaller fragments in such collisions. Colli- 
sion strength is likely to be very small for objects in the 
outer Solar System which are thought to be composed 
primarily of ices. Comets are presumably the closest 
existing analogs of distant planetesimals and they are 
known to have small internal strength, e.g. from obser- 
vations of tidal disruption of Shoemaker-Levy comet by 
Jupiter (Greenberg et al. 1995). Thus, it would be nat- 
ural to expect that collisional fragmentation triggered 
by the dynamical excitation of planetesimals by massive 
protoplanetary cores readily occurs at least in the outer 
Solar System. 

Efficiency of fragmentation is set by the collision 



timescale ti„ e i of planetesimals with sizes in which most 
of the solid mass is concentrated. Using H73JI for 1 km 
bodies with \ = 1 and v z ~ Vh we estimate 9 it to be 
~ 10 3 yr at 1 AU and ~ 10 7 yr at 30 AU. This might 
seem like a rather long timescale in the outer Solar Sys- 
tem but one should keep in mind that channeling just 
10% of mass into the population of cold planetesimals 
would increase the embryo's accretion rate by a factor of 
« 10 compared to the accretion of dispersion-dominated 
bodies, and this can be accomplished in a time 10 times 
shorter than ti ne i, i.e. in about several Myrs. Note that 
according to (|84|l the characteristic growth time of 10 25 
g embryo by accretion of small bodies is also several 
Myrs for \ w 0.1. Growth time decreases as \ grad- 
ually goes up meaning that several Myrs is a natural 
evolution timescale for such embryos in the outer Solar 
System. We may conclude that if embryos grow mainly 
by accretion of small planetesimals then the planet for- 
mation timescale is intimately related to the timescale of 
fragmentation of massive planetesimals in catastrophic 
collisions. 

Accretion rate and dynamics of protoplanetary cores 
would in the end depend on the details of the time evo- 
lution of the mass fraction in small planetesimals x> see 
^3 3.6l andlBl Scaling of x with time would also determine 
whether inelastic collisions between small planetesimals 
are an important dynamical factor; as we demonstrated 
in SJSJ this can in some cases be an issue in the inner parts 
of the protoplanetary nebula. The amount of solid ma- 
terial contained in small planetesimals is determined by 
(1) the input of mass due to the fragmentation of large 
planetesimals, (2) the removal of mass via the accretion 
by embryos, and (3) the evolution of the surface den- 
sity of small planetesimals due to the random scattering 
by embryos and their inward migration induced by the 
gas drag (Weidenschilling 1977; Thommes et al. 2003). 
Self-consistent calculation of x has to combine all these 
contributions and is beyond the scope of this study. 

Rapid accretion of small planetesimals can proceed 
provided that not only small bodies but also embryos 
themselves are on almost circular and uninclined orbits. 
As we have demonstrated in H3 3.61 oligarchic growth 
allows simultaneous existence of many protoplanetary 
cores only if cores are lighter than about 10 25 — 10 27 
g (at 1 AU), since in that case their random motions 
can be kept small by planetesimal dynamical friction. 
The exact value of the maximum core mass at which 
dynamical stability is still possible sensitively depends 
on the amount of mass concentrated in small shear- 
dominated bodies, see H3 3.BI After reaching this mass 
embryos would be dynamically excited (e, i > Rr/o) 
and even if small planetesimals can still be kept confined 
to a cold, thin disk, embryo's accretion would proceed in 
the dispersion-dominated regime rather slowly (because 
relative embryo-planetesimal velocity is increased above 
QRh). 

Discrete nature of the shear-dominated scattering of 
small planetesimals by the embryos is very important for 
accurate calculation of processes characterized by the en- 
ergy threshold such as the disruption of small planetes- 

9 Collision timescale would be smaller for small v z /vh (see Jit , 
but this is partly compensated by the burst-like character of the 
fragmentation process, see below. 
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imals in catastrophic collisions. Usual continuous ap- 
proximation suitable in the dispersion-dominated regime 
would not work in such cases because it characterizes 
planetesimal velocity by its average value which is either 
above or below the threshold, meaning that correspond- 
ing process is either always on or always off. In reality, 
planetesimal velocity changes continuously between the 
encounters with embryos from very large values (about 
QRh) to very small ones. As a result, planetesimal veloc- 
ity can be above the threshold for some time and during 
this period corresponding process would operate (see Fig- 
ure EJl- Later on velocity would drop below the thresh- 
old and process would switch off. This is qualitatively 
different from what one would obtain using continuous 
description for the planetesimal scattering in the shear- 
dominated regime. For the fragmentation of small bod- 
ies this implies that planetesimal destruction in catas- 
trophic collisions proceeds in bursts right after each pas- 
sage of the embryo, when relative planetesimal velocities 
are high, but later on, when velocities are damped by 
the gas drag, collisions might not be energetic enough to 
continue fragmenting planetesimals. 

8. CONCLUSIONS. 

We have explored the details of planetesimal dynam- 
ics near protoplanetary embryos in the presence of gas 
drag. We showed that large (> 1 km in size) planetes- 
imals are kept in the dispersion-dominated regime as a 
result of scattering by protoplanetary cores scattering, 
although their average velocities are reduced by the gas 
drag. Bodies smaller than roughly 0.1 — 1 km (depending 
on the location in the nebula) interact with protoplan- 
etary cores in the shear-dominated regime. Owing to 
the action of the gas drag these planetesimals settle into 
a geometrically thin layer near the nebular midplane; 
between consecutive encounters with the embryos they 
experience strong velocity damping which allows them 
to approach embryos every time with the relative veloc- 
ity comparable to the Hill velocity — minimum velocity 
which can be achieved in a differentially rotating disk. 

For different locations in the proto-Solar nebula we 
have determined which planetesimals are only weakly 
affected by the gas drag and are in the dispersion- 
dominated regime, and which planetesimals are so 
strongly coupled to the gas that their velocities are below 
the Hill velocity of the protoplanetary cores. Dynamical 
peculiarities of the shear-dominated regime in the pres- 
ence of gas lead to a very high efficiency of accretion of 
small bodies by the embryos. If the surface mass density 
of small bodies is high enough (> 1% of the total sur- 
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face density in solids) their accretion would dominate the 
embryo's growth rate (relative to the accretion of more 
massive, dispersion-dominated planetesimals). 

Large planetesimals (> 1 km in size) likely containing 
most of the mass in solids have large random velocities 
and are not very efficiently accreted by embryos. How- 
ever, they inelastically collide with each other at high 
velocities and fragment into smaller pieces contributing 
to the population of small bodies. Thus, the embryo's 
growth by accretion of small shear-dominated planetes- 
imals can be regulated by the fragmentation of bigger, 
dispersion-dominated bodies. Planetesimal fragmenta- 
tion would probably be easiest in the outer Solar System 
where colliding bodies are mostly composed of ices and 
are therefore internally weak and susceptible to easy de- 
struction. The natural timescale for the growth of 10 25 g 
protoplanetary embryos by accretion of small planetesi- 
mals turns out to be around several Myrs at 30 AU from 
the Sun. Cores of ice giants can be formed in ~ 10' 
yr after large planetesimals lose all their mass in catas- 
trophic collisions to small debris which can be effectively 
accreted by massive embryos. This scenario would work 
only if the population of protoplanetary cores formed as 
the outcome of oligarchic growth can be kept on almost 
non- inclined and circular orbits. As we demonstrated 
in H3 3.61 this is possible only if embryo masses do not 
exceed a specific limit dictated by the efficiency of plan- 
etesimal dynamical friction, which sensitively depends 
on the planetesimal mass fraction \ locked up in small 
planetesimals and velocity of large dispersion-dominated 
planetesimals. 

Future work should address the issues of the self- 
consistent evolution of the mass fraction \ contained in 
small bodies; the role of inelastic collisions using im- 
proved fragmentation physics (can be an issue in the 
inner Solar System); the final fate of massive embryos 
which cannot be kept dynamically cold by the planetes- 
imal dynamical friction, and so on. 
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Table 1. Core accretion and dynamical friction regimes. 



Region 10 a Conditions Dyn. friction 11 ^ Accretion 12 b M e}Cr /Mi so ls c 
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SD 18 b 


SD 
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II 
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IV 
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x 3 ? 2 (v/nR H f 



Wetherill, G. W. & Stewart, G. R. 1989, Icarus, 74, 542 Zappala, V., Cellino, A., Dell'oro, A., Migliorini, F., & Paolicchi, 

Weidenschilling, S. J. 1977, MNRAS, 180, 57 P. 1996, Icarus, 1996, 124, 156 

Whipple, S. J. 1971, in From Plasma to Planet, Almquist & 
Wiksell, p. 211 

Roman numerals correspond to different regions in Figure Shows which planetesimal population — 
dispersion-dominated (DD) or shear-dominated (SD) — dominates dynamical friction and accretion rate of embryos. 
M GjCr and M iso are given by l(55j) and (|5Ujl in regions I & II and (|5?|) and (|5l"1) in regions III & IV. 
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